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Abstract 

We study the photon-quark-quark and Higgs-gluon-gluon form factors for on-shell massless quarks 
and gluons in perturbative QCD. Previous third-order results for the quark case are extended by 
calculating the fermion-loop contributions up to the finite terms in dimensional regularization. For 
the gluon case the complete set of infrared poles at three loops is derived. Using the exponenti- 
ation of the form factor, the latter results are employed to extract a function entering the infrared 
factorization of general third-order amplitudes. We evaluate the infrared finite absolute ratio of 
the time-like and space-like gluon form factors up to the fourth order in the strong coupling con- 
stant. The result supports previous indications that the perturbative expansion of the Higgs boson 
production rate at the LHC is under control. 



The form factors of quarks and gluons, i.e., the QCD corrections to the qqX and ggX vertices with a 
colour- neutral particle X, are quantities of considerable phenomenological and theoretical interest. 
These three-point amplitudes represent important gauge invariant (if infrared divergent) parts of 
the perturbative corrections to processes of utmost relevance, like the production of lepton pairs 
and the Higgs boson H at proton colliders. Due to the exponentiation of the form factors [1-3], 
already the dimensionally regulated pole terms at order a", combined with appropriate lower- 
order information, suffice for deriving the finite absolute ratios of the time-like and space-like 
form factors to order a" +1 , thus providing valuable information about the size of the higher-order 
corrections. The form factors are also the simplest amplitudes to which the infrared factorization 
formulae of Refs. [4, 5] can be applied. Consequently explicit calculations of these quantities 
lead, through the resulting determinations of otherwise unspecified functions in these formulae, to 
unambiguous predictions for the pole structure of higher-order amplitudes involving more partons. 

In a recent article [6] we have extracted all three-loop pole terms of the electromagnetic form 
factor J ? (oc s ,<2 2 ) of on-shell massless quarks from the calculation of the third-order coefficient 
function in deep-inelastic scattering [7] and extended the resummation of this form factor to the 
next-to-next-to-leading contributions. In this letter we present the corresponding results for the 
Hgg gluon form factor f 8 ((X Sl Q 2 ). As a first step towards the complete three-loop computations, 
we moreover extend the fermion-loop (ny) parts of J q to the order £° in dimensional regularization 
with£> = 4-2£. 

Since results on deep-inelastic scattering are our starting point, we first address the space-like 
form factors. Thus the relevant amplitude for the y*qq quark case is 

= ie q (wy^)^(a s ,2 2 ), (1) 

where e q represents the quark charge and Q 2 the virtuality of the photon. The corresponding Hgg 
vertex is an effective interaction in the limit of a heavy top quark, 

L cS = ~C H HG^G a ^. (2) 

Here G^ v denotes the gluon field strength tensor, and the prefactor Ch is determined by the top- 
quark loop including all QCD corrections. Neither these corrections nor the renormalization of 
G" N G a ^ v are relevant to the present study, hence we refer the reader to Refs. [8,9] for details. 

Besides the electromagnetic case utilized in Ref. [6], we have also calculated the third-order 
corrections to deep-inelastic scattering by exchange of a scalar interacting with gluons according 
to Eq. 0. The inclusion of this (experimentally entirely irrelevant) process was required for 
obtaining the full set of splitting functions governing the evolution of the parton distributions at 
the next-to-next-to-leading order [10, 11]. All terms up to e° were consistently kept also in this 
part of the calculation. Hence we can repeat the procedure discussed in Ref. [6] for ^F 8 (a s ,Q 2 ), 
and derive the complete set of pole terms, £~ 6 . . . e , at three loops. 

In order to access the finite (e°) parts of the form factors in this approach, the three-loop com- 
putations of deep-inelastic scattering need to be extended to order £. This is not an option for the 
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complete calculation with its more than 100000 tabulated integrals. However, only very few gen- 
uine three-loop diagrams remain if we confine ourselves to the ny contributions to the quark form 
factor (these diagrams can be found in Ref. [12]). As neither of those diagrams is of a particularly 
difficult type, we were able to carry out the required extensions for this subset. 

We present our results in terms of the expansion coefficients jf of the bare (unrenormalized) 
form factors (p = q, g) 

^(a s b , 2 2 ) = l + £ (flf) f^j tf, «s-|- (3) 



The fermionic contributions to the quark form factor are given by 

7l f = 0, (4) 
a ^ f 1 14 1 /353 1 \ 7541 14„ 26,, 



/ 150125 353. 364 41 2 
+ %^944" + ^ 2 - ^ 3 - 30^ 



2 /2877653 26. 7541. 4589. 287. 2 242. , . 

+£ Uiw~t^ 3+ ^ 2 ~^ 3 ~^ 2 ~75-^' (5) 



^2 f 2 37 1 / 545 1,\ 1 / 6499 133. 

= CFn f\-3&-9& + #\- 27""3 C2 J + ^r^T"l8" C2 
146. \ 1 / 138865 2849. 2557. 337. 2 
+ ^ 3 J + £l"^86-"^ 2 + ^ 3 + ^ 2 



2732173 45235 „ 51005. 8149. 2 343 „ „ 278 



2916 162 " 81 " 216 " z 9 45 
^ ^ f 88 1 /2254 16. \ 1 /13679 316. 

». \ 1 / 623987 11027. 6436. 44. 2 \ 
<3 ) + i (,-2187" + "W^ 2 ~ "gT^ 3 ~ T^ 2 ) + 



27 7 e \ 2187 243 ' 81 ' 5 ^ L ) 6561 
442961 45074 1093. 2 368. 208. 

+ ~H5^ 2 — sT * _ "IrT^ 2 + "^ 3 " T^ 5 

2 / 8 188 1 / 124 4 \ 1 / 49900 94 
+ C F n f\ "8T?~ 243? + ?V '27"~9 fa J + iV "W ~ 27 = 



136. \ 677716 62. 3196. 83 y2 \ 
+ Jl^ 3 ) ~ "656T " T^ 2 + W^ 3 " 135 C * J • (6) 

Here stands for the number of effectively massless quark flavours, Cf and Ca are the usual QCD 
colour factors, Cf = 4/3 and Ca = 3, and the values of Riemann's zeta function are denoted by C, n . 
The E 1 and e 2 terms in Eq. © and the corresponding non-fermionic contributions [6] have recently 
been confirmed in Ref. [13], where exact expression for both two-loop form factors were derived. 
The finite (e ) contributions in Eq. © represent a new result of the present letter. 
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The first three expansion coefficients © of the unrenormalized gluon form factor read 



7. , 62„\ A { 14 949 3 49 2 
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-y^3+32^ 5 
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481 209 247 2 \ 1 / 145219 20329 241 3751 2 

+ ^ 2 "^ 3 + ^ 2 J + 8r^T^ + ^4T^ 2 + ^ 3 "W^ 2 

85,, 878, \1 . f 2 2 1 /1534 47, \ 
- y ^3 - I + Q.^ - ^ - ^ + ? + -C2J 

1 /4280_ 425 518 \ 1 /_ 92449 _ 7561 1022 2453 2 
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[ 20 1 /526 160, \ 1 /2783 22, 224, 176, 2 \ 1 
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+ C ^\3i| +C ^\-^- 243? + ?U + ^ 2 J + iVW + ^ 2 

+ ^3j}+C F n/|^ + -^--yC 3 j}. (8) 

As in Eq. ©, the one- and two-loop quantities Q and © are written down to the accuracy in £ 
required for the extraction of the three-loop form factors to order 6 . The terms up this order in 
Eq. dU) have been obtained before in Refs. [8,9]. Our corresponding coefficients of E 1 and £ 2 agree 
with the recent all-order expression of Ref. [13]. Eq. ®, of which we will present an independent 
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check below, is the main new result of this article. 

In the context of Refs. [14, 15] it is interesting to note that the terms of highest transcendentality, 
i.e, the coefficients of C, n and C,iC,j with i + j = n, in the £~ 2 '+" contributions to Ji agree between 
the quark and gluon form factors for the Super- Yang-Mills case Ca = Cp = n c . In fact, up to two 
loops this can be readily shown to all orders in £ by expanding the prefactors of the master integrals 
inEqs. (8) - (11) of Ref. [13]. 

The exponentiation of the form factors is performed for the (coupling-constant) renormalized 
quantities. These are obtained from Eqs. © - © by replacing the bare coupling a* by its renor- 
malized counterpart a s according to 

*2 i R \ /R3 



where Bi are the usual expansion coefficients of the beta function of QCD, Bo = -y C4 — §«y etc. 
Note that, unlike Refs. [9, 13], we do not include the multiplicative renormalization of G£ v G a ' pV in 
Eq. © into the definition of the (renormalized) gluon form factor. 

The exponentiation of the (coupling-constant) renormalized form factors can be written as [2,3] 

2 e 2 /<" 2 1 

ln^a s ,^,£^) = X - J ^\K(a s ,E)+G(l,d(^a s ,E),E)+J ^A(a(X,a s ,E))J . (10) 

Here K p (a s ,e) are scale-independent counter-term functions consisting of a series of poles in £ 
in the MS scheme. These functions can be determined recursively from a renormalization group 
equation (see, e.g., Ref. [3]) in terms of the cusp anomalous dimensions A p [16]. The functions 
G p , on the other hand, can be expanded in non-negative powers of £ at all orders of a s . Finally a 
is the running coupling in D dimensions, see Refs. [6, 17]. After expansion in a s , the integrals in 
Eq. (TTOb can be solved using algorithms for the evaluation of nested sums [18, 19], see Ref. [6]. 
Transforming back to the unrenormalized expansion coefficients ff in Eq. © the results read 



7i = ^A? + ilAi(2Gi-Po) + il(^-A 2 -2PoGi)-~G2, (12) 



^3 = -^^A^11a?(G 1 -Bo)- t ^1a 1 (9G?-9A 2 -27B G 1 + 8B2) 

1 l -{3G\-9A 2 Gi - I8A1G2 + 4Mi " 18B Gf + I6M2 + 24B^Gi; 



144 £ 3 

+ 72? (9GlG2 " 4A3 ~ 6plGl ~ 24poG2) ~ \\ G * ■ (13) 
The corresponding expression for the four-loop form factors can be found in Ref. [6]. 

From now on we confine ourselves to the gluon form factor ? g . The corresponding anomalous 
dimensions A 8 are related by a factor C A /C F to the quark quantities A q [16] which are known to 
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order a s 3 [ 1 0, 20] . This relation has been verified by the explicit calculation of Ref . [11]. As always 
using the expansion parameter a s = oc s / (47l), the available coefficients read 



Af 


= 4C A , 


Af 


= 8Ci( 


Af 


= 16Cj 



(14) 

^Z-^+SQn/-^, (15) 



245 67 . 11 . 11 . 2 \ ^ ^ / 55 „ „ 

+ 16C ^/ (-M + y C2 " I ^ 3 ) + l6C ^ (-27 ) • (16) 

Beyond the third order only the (small) leading-n^ contributions are known [21]. 

Using these coefficients, the /-loop pole terms £T 21 . . . e . in Eqs. dTTT) - <fT5T> can be predicted 
from lower-order information, thus providing a strong check of either the exponentiation formula 
or the explicit /-th order calculation. Our results pass this check, thus Eqs. Q - © can be em- 
ployed to recursively derive the resummation functions Gf for i = 1, 2, and 3 to the respective 
fifth, third and zeroth order in e, yielding 

Gf = £C A (-2^ 2 ) +e 2 C A (4-^+£ 3 C A (l2-^2 2 )+e 4 Q(28-2^ 
14. . 124. \ . / _ 28 r 949 r3 , 98 r2 



+ -^3 - — Cs ) + e 3 Q ( 60 - 6^2 - yCa - ^C! + y £ ) , (IV) 



9022 134 „ 88„\ „ /3448 20^ 80>, 



81 3 3 3 'V A J V 81 3 3 3 



^ /^ 134 2^2/^ 141 677 568^ 4556^ 671„ 2 
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V 54 9 

^ „i( 373975 27320. 4096. 1276. 2 80. . 
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„, /266072 7328,. 56, 328, ,\ „ „ Z3833 „, 



752 „ 32 
"5 



^ ,^2 ^ 2/^ 28114 160 >- 256. \ 

-4C F n f + C A nf[ 1 2 1 3 

J \ 72 9 27 Si 27 S3 / 



+ C W (_104 + « 6 ). (19) 

Results analogous to the terms up to order £ in Eqs. ( fT7t and ( fT8l have been derived before, in a 
different notation, in Ref. [9]. The third-order coefficient (fT9l) is a new result based on Eq. ©. 
Note that G| is one of the quantities which determine the infrared structure of QCD amplitudes in 
the framework of Refs. [4, 5] . 

We now turn to the check of our main results © and (IT9"1) announced above. Inspired by a key 
observation of Ref. [9], we decompose the resummation functions Gf (p = q, g) according to 

Gf = 2(*f-8 M Po) +/f° + eGf , 

Gf = 2(5 2 p -28 pg pi)+// + (3 Gf(e = 0)+£G 2 D , (20) 
Gf = 2(fif-38 pg p 2 )+/ 3 p + p 1 Gf(e=0) + Po[Gf(8=0)-PoGf(6=0)]+6Gf 

with 

F = z~ l [F-F(e=0)] (21) 

and 5,f denoting the coefficients of 8( 1 —x) in the n-loop diagonal MS splitting functions Pp£~^ (x) 
[10, 1 1]. The crucial point of this decomposition is that the functions ff are, like the cusp anoma- 
lous dimensions A„, universal up to the factor C A /C F , i.e., ff = C A /C F ff. The functions ff also 
exhibit the same maximally non-Abelian colour structure as the A„ [16], with 



ft 


= o, 


f! 


= ^cJ 




= c F c\ 



(22) 



56 „ /404 



-PoC 2 --n / + C 4 (— - 14^ ) | , (23) 
/136781 12650^ 1316 „ 352 ^ 2 176 „ „ inn „ 



11842 2828 „ 728 „ 96 „,\ , / 1711 
304 „ 32„,\ „ ,/ 2080 40 „ 112 



+ C A C F n f ( 1 C 2 H C3 to +C F n fl 

f \ 729 81 27 S 5 S2 / F A 27 



+4fc + ^?3 + f gj+C^-^-^fc+^j . (24) 

Consequently, once A 9 is known from the explicit calculation of the three-loop quark form fac- 
tor [6], the pole terms of the gluon form factor can be derived from the results of Refs. [10, 11] 
and lower-order quantities. This procedure confirms Eq. (IT9l and thus Eq. ©. For applications 
to other cases, like the coupling of a preudoscalar Higgs boson, it should be noted that the terms 
subtracted from Bf in Eqs. (l20l are given by the renormalization of the operator in Eq. ©. 
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As a first concrete application of our results we consider the absolute ratio 1 7 8 (q 2 ) / T 8 {—q 2 ) \ 
of the renormalized time-like and space-like form factors. This quantity is infrared finite and 
directly enters the cross section for Higgs boson production in hadronic collisions. In terms of the 
coefficients A, and G, (£ = 0) this ratio is given by [6] 

l+a s {3^ 2 A l } + a 2 {^A 2 + 3^oG l +A 2 )} 
+ a s 3 \^QA\ + 3^ 2 A! (3p Gi - p 2 + 3A 2 ) + 3^(A 3 + pi G, + 2f3 G 2 ) } 

+ 4{ f C 2 X + ^ 2 3 A?(3PoGi - 2f3 2 + 3 A 2 ) + |g(-6PjjA 2 + 3f3 2 G 2 
+ 3A\ + 12(3oAiG 2 + 6A1A3 + 6MiGi " 5poPiAi + 6P0A2G1 
-6p^Gi)+3C 2 (A 4 + |32Gi+3|3oG3 + 2piG2)} + o(a s 5 ) (25) 

for the couplings a s (q 2 ) = a s (—q 2 ) = a s . All terms contributing at the fourth order are now known, 
with the exception of the small four-loop cusp anomalous dimension A|. For the latter we can 
employ the [1/1] Pade estimate of the quark case in Ref. [22], accordingly multiplied by C A /C F , 

A\ w 17660, 9704, 3949 for ^ = 3,4,5. (26) 

Below a conservative (but numerically irrelevant) uncertainty of 50% is assigned to this estimate. 
Inserting the numerical values of Eqs. (fl4l - (fT9l for = 5 quark flavours, we obtain 

2 

= l+4.712a s + 13.69a s 2 + 25.94a s 3 + (36.65±0.35)a s 4 . (27) 

The contributions up to order a 2 have been derived before in Ref. [8], where the factor 153.2 in the 
penultimate line of Eq. (19) should read 135.2. The low-order contributions to Eq. (l27l) are larger 
than for the corresponding quark quantity in Ref. [6]. On the other hand the growth of the coeffi- 
cients is slower in the present case, surprisingly suggesting a more benign higher-order behaviour. 
As for Higgs boson production couplings of the order oc s ~ 0. 1 are relevant, Eq. (l27l) further sup- 
ports indications from the next-to-next-to-leading logarithmic threshold resummation [23,24] that 
the corrections beyond the next-to-next-to-leading order (a 2 ) of Refs. [25-28] are rather small. 
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